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ABSTRACT 



Aims. A strong, hard X-ray flare was discovered (IGR J12580-I-0134) by INTEGRAL in 2011, and is associated to NGC 4845, 

a Seyfert 2 galaxy never detected at high-energy previously. To understand what happened we observed this event in the X-ray band 

on several occasions. 

Methods. Follow-up observations with XMM-Newton, Swift, and MAXI are presented together with the INTEGRAL data. Long and 

short term variability are analysed and the event wide band spectral shape modeled. 

Results. The spectrum of the source can be described with an absorbed (A'h ~ 7 x 10"" cm"") power law (T == 2.2), characteristic 

of an accreting source, plus a soft X-ray excess, likely to be of diffuse nature. The hard X-ray flux increased to maximum in a few 

weeks and decreased over a year, with the evolution expected for a tidal disruption event. The fast variations observed near the flare 

maximum allowed us to estimate the mass of the central black hole in NGC 4845 as ~ 3 x 10^ Mq. The observed flare corresponds to 

the disruption of about 10% of an object with a mass of 14-30 Jupiter. The hard X-ray emission should come from a corona forming 

around the accretion flow close to the black hole. This is the first tidal event where such a corona has been observed. 

Key words. X-rays: galaxies - X-rays: individuals: IGR J12580-I-0134 X-rays: individuals: NGC 4845 



1. Introduction 

NGC 4845 is a nearby, high surface brightness spiral galaxy 
with morphological type SA (s)ab sp classified as a Seyfert 2 
dVeron-Cettv & Veronll2006h . The distance to NGC 4845 was 
determined as Z) - 15 . 6 Mpc based on the TuUy-Fisher relation 
(iTullv & Fisheilll988h . lCorsini et al.l (Il999h obtained D = 13.1 
Mpc by including the influence of a dark matter to the systematic 
velocity. We adopt in this paper D = 14.5 Mpc {Hq = 75 km s ' 
Mpc ') after IShaplev et al] (1200 ih who performed a thorough 
literature search to determine the most reliable distance of this 
galaxy. 

The galactic disk inclination is 76° dPizzella et alJ 120051) . 
Coordinates of the galax y (i.e. g.i7nnn - 12' ' 58'"0 1"19 and 
feooo = +01°34'33'.'02. IVeron-Cettv & VeronlF2010l) set the 
Galactic neutral hydrogen column d ensity into the direction of 
the source to Ng ' = 1.67xl0^"cm ' dPickev & Lockmanll990l: 
IStark et al.|[T9^ . 

NGC 4845 has been observed, so far, in the optical 
(e.g. by the Sloan Digital Sky Survey - SDSS, the Palomar 
Optical Sky Survey - POSS and the Hubble Space Telescope 
- HST), and in the infrared domain (by the Two Micron 
All Sky Survey - 2MASS and the Infrared Astronomical 



Satellite - IRAS) JM pshir et a ll 119901: ISp inoglio et al 
Diorgovski et al. 19985 ISanders et al.ll2003[ iJarrett et all 
Schneider et al...2005.) . 



1995; 



2003; 



The galaxy has also been observed by the Green Bank 
Telescope and the Very La rge Array in the radio bands 1.4, 
4.8, 8.4 GHz (iCondon et al., 1998, 2002) The radio emiss ion is 
dominated by star-forming regions (see iFilho et all 2000[ their 
Fig. 5d) and does not reveal water maser emission (iBraatz et alJ 



l2003h . Monitoring in the X-ray domain was carried out using 
the Imaging Proportional Counter onboard the EINSTEIN satel- 
lite dFabbiano et al.|[l992l) . The authors estimated only an upper 
Umit in the 0.2-4.0 keV energy band of Fx < 2.52 x 10 '^ erg 

-2 -1 

cm s . 

INTEGRAL disc overed a new hard X-ray source 
IGR J12580"H0134 (iWalteret al.ll201lh during an observation 
performed in the period January 2-11, 201 1, with a position (i.e. 
RA = 194.5212 deg DEC = 1.5738 deg, +2.3 arcmin, J2000.0) 
consistent with that of NGC 4845. A few days later, Swift/XRT 
and XMM-Newton observations confirmed the association with 
the central regions of the Seyfert 2 galaxy. IGR J12580-H0134 
was detected at a peak flux F2-iokeV > 5.0 x 10"' ' erg cm"' s"', 
corresponding to a variability by a factor > 100, which is very 
unusual for a Seyfert 2 galaxy. 

In sect. |2] we report the analysis of the available high- 
energy observations of IGR J12580-H0134 performed with 
XMM-Newton, Swift, MAXI, and INTEGRAL. We analyse these 
results in sect.[3]and discuss them in sect. |4] 



2. High-energy observations of IGR J1 2580+01 34 

We analysed INTEGRAL, MAXI, Swift, and XMM-Newton obser- 
vations of IGR J12580-I-0134. A log of the pointed observations 
is given in Tables[T]and|2l Throughout this section, all uncertain- 
ties are given at 1 sigma confidence level. 
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Table 2. INTEGRAL observations of IGR J12580+0134. 



Revolutions 


Start time 




Stop time 


EXP 




(UTC) 


(MJD) 


(UTC) 


(ks) 


1004-1006 


2011-01-02 13:51:14 


55563.57725 


2011-01-11 04:47:40 


100.0 


1009 


2011-01-18 19:09:42 


55579.79841 


2011-01-22 01:49:00 


1.6 


1012 


2011-01-27 10:55:12 


55588.45500 


2011-01-27 11:21:46 


1.0 


1052-1055 


2011-05-28 00:16:37 


55709.01155 


2011-06-06 03:15:17 


44.2 


1057-1061 


2011-06-10 03:14:24 


55722.13500 


2011-06-24 03:55:16 


90.4 


1063 


2011-06-28 12:02:32 


55740.50176 


2011-06-30 17:10:01 


29.0 


1067-1068 


2011-07-10 21:03:15 


55752.87726 


2011-07-15 16:05:07 


45.0 


1070 


2011-07-20 08:31:40 


55762.35533 


2011-07-21 16:29:40 


11.7 



Notes. EXP indicates the effective exposure time. 



Table 1. Swift and XMM-Newton observations of IGR 
J12580H-0134. 



Instrument 


Obs. ID 


Start time 


Net exp. 






(UTC) 


(ks) 


Swift 


00031911001 


2011-01-12 23:59:01 


3.0 


Swift 


00031911002 


2011-01-13 03:26:01 


2.1 


XMM-Newton 


0658400601 


2011-01-22 16:23:28 


14.0 


Swift 


00031911003 


2012-06-29 07:56:59 


3.3 



2.1. INTEGRAL 

INTEGRAL data were analysed using OSA software (ver 9). We 
considered data from IBIS/ISGRI in the 17.3-80 keV. The aver- 
age fluxes and spectrum were extracted from the mosaic images 
with mosaic_spec. 

INTEGRAL observations are divided into pointings with du- 
ration of ~ 2-3 ks, which are collected during the satellite revo- 
lutions. We considered in this paper all available (i.e. 555) point- 
ings from the discovery observation to July 21, 2011, included 
in revolutions 1004-1070. 

Figure [H show s the IBIS/ISGRI ( iLebrun et alj l200l 



ire m 
let alii 



lUbertini et ani2003h mosaic obtained around IGR J12580-H0134 
(17.3-80 keV, significance map) accumulated from January 2 to 
11, 201 1. During this period the source significance reached 9.7 
sigma in the 20-40 keV energy band, corresponding to a flux of 
■P'20-40keV - (2.5 ±0.3) X 10"" erg cm"" s"' for an ex posure time 
of 40 ks. The source was also detected by JEM-X ( iLund et all 
l2003h (3-20 keV) for an effective exposure time of 13 ks with a 
flux of F3-iokeV = (5.0 + 1.4) X 10-" erg cm^^ s"'. 

We also analysed all publicly available IBIS/ISGRI data ob- 
tained on the field from the beginning of the m ission up to 
January 27, 2011 using the HEAVENS interface (IWalter et al.l 
120101) with OSA9. The INTEGRAL observations performed in 
previous years (i.e. from 2003 to 2010) did not reveal any hard 
X-ray activity, and we obtained a five-sigma upper limit of 
5.9 X 10"'" erg cm"^ s"' in the 20-40 keV energy band. 

The resulting count rates are listed in Table |3] and plotted in 
Fig.m 

2.2. MAXI 

NGC 4845 has been observed by the Monitor of AU-sky X-ray 
Image (MAXI), attached to the Japanese Experiment Module on- 
board the International Space Station. Among the two X-ray in- 
struments onboard MAXI, we analysed only the Gas Slit Camera 
i GSC) data in the present paper, because o f its higher sensitivity 
Matsuoka et al.ll2009HMihiu-a et alJl201 ll) . 



■2r.-Q- 



19^.0 




^1-5— 



194.5 



194.0 



Fig.l. The INTEGRAL IBIS/ISGRI mosaic of 
IGR J12580+0134 (17.3-80 keV, significance map) observed 
2-11 January 2011. The dash black grid denotes equatorial 
J2000.0 coordinates in degrees. 

Table 3. Intensities of IGR J12580+0134 observed by 
INTEGRAL/lSGRl. 



START TIME 


END TIME 


Count rate 


Signific. 


(MJD) 


(MJD) 


(count s"') 




During flare (i.e. Jan 2011) 






55542.63936 


" 55563.22030 


0.192 ±0.120 


1.6 


55563.57725 


55566.23899 


0.635 ±0.139 


4.6 


55566.56943 


55569.20899 


1.132±0.153 


74 


55569.56013 


55572.19978 


1.157 ±0.167 


6.9 


55579.79841 


" 55588.47273 


1.652 ±0.674 


24 


During flare - 


total 






55563.57725 


55572.19978 


0.944 ± 0.087 


10.8 


After flare (i.e 


. Jun-Jul 2011) 






55713.98265 


55723.36030 


0.341 ±0.125 


2.7 


55725.96026 


55742.71529 


0.304 ± 0.083 


3.7 


55752.87726 


55763.68727 


0.013±0.111 


0.1 


After flare - total 






55713.98265 


55763.68727 


0.269 ± 0.056 


4.8 


Period < 2011 (Jan 2003-Dec 2010) 






52644.50424 


" 55545.30596 


0.014 ± 0.026 


0.5 


Period > 2011 (Jan-Jul 2011) 






55563.57724 


55763.68724 


0.470 ± 0.047 


9.9 



Notes. " Data obtained through the HEAVENS interface. 
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Fig. 3. The optical finding chart of NGC 4845 taken from 
the Space Telescope Science Institute - Digitized Sky Survey 
(STSci-DSS). White soHd contours show the XMM-Newton 
Epic-pn image of IGR J12580+0134 (0.8-10 keV, significance 
map). The black grid denotes equatorial J2000.0 coordinates in 
degrees. 



Fig. 2. The light curve of IGR J12580+0134 observed in the 
17.3-80 keV energy band by INTEGRALflSGRl (red points) 
and in the 2-20 keV band by MAXI (blue points). The ar- 
row shows an upper limit. Red short-dash line indicates the 
5cr upper limit of the source observed by IBIS/ISGRI before 
2011. Red long-dash line shows an upper limit observed by 
EINSTEIN ( lFai7bianoetal.lll992h and interpolated to 17.3-80 
keV. The interpolation is based on the model fitted to the XMM- 
Newton/ INTEGRAL observation. The interpolated EINSTEIN 
flux is very likely overestimated because of the harder spectrum 
observed in 2011. 



visible sources area, was also chosen as a circle with a radius of 
50 pixels. There are a number of bad columns on the XRT CCD, 
and the corresponding pixels were not used to collect the data. 
To account for them, we used a corrected exposure map created 
with the XRTEXPOMAP command and used it to create the ancil- 
lary response files with the command xrtmkarf . We maximised 
the signal-to-noise (S/N) by summing up all the available data 
obtained in January 2011 (effective exposure time 5116 s). The 
final extracted spectrum was rebinned to have a minimum of 30 
counts in each energy bin. Only 12 source photons were detected 
in June 2012. 



Figure |2] shows the MAXI Ught curve of NGC 4845 in the 2- 
20 keV band, obtained from the MAX//GSC on-demand process 
interface with the source and background radii of 2 and 3 deg, 
respectively, centred on the galaxy. 

2.3. Swift 

The field of view (FOV) around IGR J12580+0134 was ob- 
served three times with Swift/XRT on January 12 and 13, 
2011 and on June 29, 2012 for a total exposure of 3.1, 
2.1, 3.3 ks (observations ID 00031911001, 00031911002, and 
00031911003, PI Walter), respectively. We analysed data col- 
lected in photon-counting mo de (PC) by using standard pro- 
cedures dBurrows et al.l l2005b and the latest calibration files. 
The Swift/XRT 'enhanced' position of the source is ffjaooo - 
194.50492 deg and 5j2ooo = 1.57579 deg with the error 
of 2.1 arcsec. To o btain those values we followed the pro- 
cedure described in lEvans et all (|2009|) where we used the 
Swift/XRT-Swift/VWOT aUgnment and matching UVOT FOV to 
the USNO-Bl catalogue. We extracted light curves and spectra 
of IGR J12580H-0134 using die Level 2 event filefl No pile-up 
problems were found in either observation. Therefore, we used 
circle region with radius of 20 pixels (~47 arcsec) around the 
source coordinates. The background region, which lies in the no 



' The online XRT 

|http://www.swift.ac.uk/XRT.shtml| 



analysis 



threads 



2.4. XMM-Newton 

XMM-Newton observed IGR J12580-H0134 on January 22, 201 1 
(obs. ID 0658400601, total exposure time ~ 21 ks, PI Walter) 
and obtained a source position of ajaooo = 194.50427 deg and 
feooo = 1.57585 deg 4.1 arcsec (Fig.O. 

During this observation all the EPIC-pn and EPIC-MOS 
cameras were operated in full frame mode. The galaxy was ob- 
served with both the EPIC-pn and EPIC-MOS; however, in this 
paper we focus on the data obtained from the first detector 

We processed XMM-Newton observation data files (ODFs) 
with pipeline eproc (Science Analysis System, SAS, v. 10.0) in 
order to produce a calibrated event list. The event file was fil- 
tered to exclude high background time intervals following the 
recommendations of the SAS analysis threacQ. We extracted the 
light curve in the 10-12 keV energy band using all field of view 
(FOV). We excluded from further analysis time intervals during 
which the count rate in the 10-12 keV energy band was higher 
than 14.0 count s"'. The resulting effective exposure time was 
14.02 ks. We extracted spectra and relevant light curves in the 
0.1-2,2-10, andO.l-lOkeV energy bands. Since the X-ray count 
rate of IGR J12580H-0I34 is high, we corrected the product for 
pile-up extracting the source in an annulus with an inner radius 
of 1.25 pixels. The background extraction region was chosen at 
~ 9.5 arcmin from NGC 4845, far away from any sources but 



How 



to 



extract 



PN 



spectra 



also 



|http://xnnm.esa.int/sas/current/documentation/threads/PN-spectrum_thread.shtnri 
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Fig. 4. XMM-Newton Epic-pn background-subtracted light 
curves of IGR J12580+0134. Top and middle panels show 
light curves extracted in the energy bands 0.1-2 keV (Soft) 
and 2-10 keV (Hard), respectively. The hardness ratio, defined 
as HR - (Hard - Soft)/(Hard + Soft), is reported in the bottom 
panel. The time bin is 100 s. 



still in the same CCD of the EPIC-pn camera. All EPIC im- 
ages and spectra were corrected for Out-of-Time (OoT) events. 
EPIC-pn spectra were rebinned before fitting so as to have at 
least 200 counts per energy bin. 



3. Data analysis and results 

3.1. Source spectrum 

We extracted the XMM-Newton Epic-pn light curves in the 0.1-2 
keV (Soft) and 2-10 keV (Hard) energy bands with time bins of 
10 s. These two light curves (rebinned to 100s) and the hardness 
ratio HR - (Hard - Soft)/(Hard + Soft) are shown in Fig. |4] 
where a short time-scale variability is clearly visible. The vari- 
ability time scale t between two points characterised by the 
highest difference in their count rates in 2-10 keV is equal to 
90 + 5 s (AAmpl= 10.5 + 2.0 count s"', found between 17 and 
19.5 ks in Fig. |5]l. The count rate in the 2-10 keV band is 6.33 
count s^', in average, with an observed standard deviation of 
2.02 count s~'. It varies between 1.08 + 0.48 and 15.33 + 1.77 
count s"'. 

The hardness ratio does not show any significant variabil- 
ity suggesting that the spectral shape does not vary with flux 
on short time scales. To check this we extracted two spectra 
by selecting time intervals in which the count rate in the 2-10 
keV energy band is < 6.3 and > 6.3 count s '. The first spec- 
trum (with the effective exposure time 5.9 ks) could be fitted by 
an absorbed power law. The resulting spectral parameters are 
Nh = (7.43 + 0.14) X lO^^cm-- and F = 2.32 + 0.04 for a 
2-10 keV X-ray flux of 4.95 x 10"" erg cm"^ s"'. The spec- 
trum extracted at a higher count rate (effective exposure time 
8.1 ks) could be described using the same model, providing 
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Fig. 5. Zoom of the XMM-Newton Epic-pn background- 
subtracted light curve showing the most variable part. The time 
bin is 10 s. 



2.13 + 0.03 for an 



erg cm ^ s ' . 



Nh = (7.11 + 0.10) X lO^-'cm-^ and F 
2-lOkeVX-ray fluxof6.94x 10" 

These spectral parameters indicate that the source spec- 
trum did not change significantly; therefore, we also extracted 
the average spectrum of the complete XMM-Newton observa- 
tion (Fig. |6l bottom) that can be represented by the parameters 
Nh = (7.21 + 0.08) x lO'^cm^^ and F = 2.19 + 0.03. The aver- 
age 2-10 keV X-ray flux is 6.09 x 10"" erg cm"^ s"'(see also 
Tablegli. 

Swift observed IGR J12580-I-0134 almost one week earlier 
(i.e. on January, 13). We performed a similar analysis to the 
one for XMM-Newton and found consistent results (Fig.|6] upper 
panel). The extracted 1-10 keV spectrum could be represented 
by Nh = (6.53 + 0.38) x lO^^cm^^ and F = 2.36 + 0. 14 for a 2-10 
keV X-ray flux of 4.97 x 10"" erg cm^^ s"'. 

The IBIS/ISGRI spectrum of the source in the 17.3-80 
keV energy bands was extracted from observations performed 
January 2-11 (eff'ective exposure of 100 ks). The 17.3-80 keV 
spectrum could be fitted (x^ by d.o.f. - 0.40) by using a power- 
law model with F - 2.17 for an average hard X-ray flux of 
F(17.3 - 80keV) = 5.44 x 10"" erg cm"^ s"'. A common fit 
of the ISGRI and XMM-Newton data (with free normalisation) 
could be obtained, leading to Nh = (7.39 + 0.10) x lO^^cm"^, 
F = 2.22 ± 0.03, and F(17.3 - 80keV) = (5.408 + 0.015) x 10"" 
erg cm"^ s"' (Fig.©. 



3.2. Intrinsic diffuse emission 

The XMM-Newton/Epic-pn spectrum shows an X-ray excess be- 
low ~ 1.2 keV (Fig. |6]l, which can be fitted by a black body 
absorbed by the Milky Way leading to a temperature of T = 
0.33+0.04 keVandafluxofFbbody(0.01 - lOkeV) = 1.12x10"'^^ 
erg cm"" s"'. This black body is not a physical model but 
does provide a good representation of the excess. Its flux is in 
good agreement with the diffuse emission radiated by galaxies 
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Table 4. The best-fit parameters to the Swift and XMM-Newton observations in January 20 11 . The fit to the 1-10 ke V energy range. 



Satellite/ 

-^ exposure 


Data 


(10^2 cm-2) 


r 


(10- 


Pl. Norm, at 1 keV 
' Photons keV"' cm"^ s"') 


^2-lOkcV 

(10"" ergcm"^ s"') 


x'IA.o.f. 


SwiftlS.Vks 


All data 


6.53 ±0.38 


2.36 ±0.14 






5.69 ± 1.33 


4.968 ±0.130 


78.8/81 


XMM/5.9ks 
XMM/8.1ks 
XMM/ 14ks 


Count rate < 6.3 
Count rate > 6.3 
All data 


7.43 ±0.14 
7.11 ±0.10 
7.21 ±0.08 


2.32 ± 0.04 
2.13 ±0.03 
2.19 ±0.03 






5.50 ±0.39 
5.54 ±0.28 
5.45 ± 0.23 


4.945 ± 0.077 
6.937 ± 0.023 
6.088 ± 0.020 


90.9/81 
193.7/150 
349.5/225 



3.3. Duration and frequency of tine fiare 




^ V t ^ t ^--*- ) + ) I - t -^ =i= 



^¥W^^^^^ 



Energy (keV) 
XMM/Epio-pn + INTEGRAL/ISGRI 




Fig. 6. Unfolded spectrum of IGR J12580-I-0134 from Swift 
observations made in PC mode (top) and unfolded XMM- 
A^e wf on/Epic -pn and /A^7:EG/ML/ISGRI spectra (bottom). The 
best-fit model (an absorbed power law in all spectra plus an un- 
absorbed black-body emission below 1 keV in case of Epic-pn 
observation) is shown by solid lines in both upper panels. Lower 
panels indicate the residuals (A;^). 




Because the source was still detected significantly in July 2011 
by INTEGRAL, the flare probably lasted for 150 days at a level 
above 2 x 10"" erg cm"^ s"'. MAXI observed the source con- 
tinuously, but at a lower sensitivity. Therefore the flare was just 
detected by MAXI close to the peak in January 201 1 . 

Such a flare of IGR J12580+0134 with a duration of 150 
days has never been detected in the past ten years by INTEGRAL, 
and it could have been missed only between August 2006 and 
December 2007 when INTEGRAL did not observe that region of 
the sky. The flare of NGC 4845 is therefore exceptional, and no 
such flare was detected for a continuous period of at least 1300 
days. 

3.4. Biack-hoie mass 

The mass of the central black hole, M^h, of NGC 4845 can be es- 
timated using methods based on the X-ray variability time scale. 

The upper limit of Mbh can first be estimated assuming that 
the shortest time variability t is related to the innermost stable 
circular orbit (ISCO). The XMM-Newton light curve observed 
in 2-10 keV band indicates t < 90 + 5 s. Additionally, assum- 
ing a non-rotating black hole, we can write that ct = Rsoiw = 
2G/Wbh/c~. The observed variability time scale corresponds to 
Mbh < 9.6 x lO^M©. 

The second technique is based on the X-ray excess vari- 
ance measurements. This method uses the relationship between 
the black-hole mass and the X-ray variability Mbh = C{T - 
2Af)/cr^xj, where T is the duration of the X-ray light curve and 
At is the bin size, both in seconds, and cr^^j, = Zf l|[(-^i ~ x)^ ~ 
'l ^A/jNxP') is the normaUsed excess varian ce (iNandra et al.l 
19971 IVaughan et al.ll200l lO'Neill et al.]|2005h . We divided the 
2-10 keV XMM-Newton light curve into two parts with length 
T~ 4 ks and ~10 ks in order to avoid the gap of 5 ks in the light 
curve. The estimated mean normalised excess variance, (cr^^J, 
of the 10 s bin light curve is 0.12;^qq^. Thus, the calculated cen- 
tral black hole mass in NGC 4845 is Mbh = 2.3![|^ x IO^Mq, 
where we assumed C = 1.42 (i.e. that the mass of Cyg X-1 is 
equal to 14.8 + 1.0 Mrr,. lOrosz et al.ll201 Jl) . The errors on the 
mass were es timated by performin g Monte Carlo simulations 
(see details in iNikolaiuk et al.ll2006h . The systematic errors, re- 
lated to the calibration of the excess variance versus mass rela- 
tionship, are larger an d the black-hole ma ss could range in the 
interval 10'* - lO^'M© (lO'Neill et al.ll2005l) . 



dBogdan & Gilfanovl201 UlJia et al.l2012h . and not very far from 
the EINSTEIN upper limit and to the residual emission detected 
by Swift XTE in June 2012 (Fbbody(0.5 - 2keV) = 8 x 10""* erg 



,-2 ^-l\ 



This soft excess is completely unrelated to the nuclear flar- 
ing activity since the nucleus is very strongly absorbed. This 
emission comes from a region outside of the dusty torus of 
NGC 4845. 



4. Discussion 

4.1. Piienomenology 

We first consider the possibility that the observed variability is 
driven by absorption, i.e. that a hole opened on the line of sight in 
the AGN absorbing to rus. The [O III] Une intensity in NGC 4845 
is 5.40 X 10"'^ erg s' dHo et alJll997h . According to the X-ray - 
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Fig. 7. Increase in x^ related to an additional soft X-ray com- 
ponent (black-body with temperature of 0.2 keV) increasing in 
luminosity, when fitted by the same model as used in figure|6] A 
soft X-ray component stronger than 5 x IC*^ erg s ' should have 
been detected, if present. 



[O III] ^5007 relationship (iPanessa et al.ll2006l) the estimated un- 
absorbed 2-10 keV flux should be ^ 6 x 10"'^ erg s"' . However, 
since the X-ray flux observed in January 201 1 was almost 1000 
times larger (Table |4}, the absorption explanation is ruled out. A 
strong increase in the X-ray flux, by a factor 10^, is required. We 
also note that the source is not Compton thick, {Nh ~ 7 x 10^^ 
cm"^, see Table |4]i, indicating that it should be bright at hard 
X-rays, when active. 

Supernova ex plosions and expanding winds may emit 
hard X-ray flares illmmler & Lewinll200l Ichandra et al.ll2009t 



ISaxton et al.l 120121) . Supernova peak X-ray luminosities have 
been observed up to lO"*" erg/s with a decline usually fol- 
lowing a f-(05-i) scajjjjg Tfhe peak X-ray luminosity of 
IGR J12580+0134 is 100 times brighter, and it declined much 
faster than observed in supernovae. This, together with the po- 
sition of IGR J12580+0134 at the very centre of NGC 4845, 
makes the supernova interpretation of the observed flare very 
unlikely. 

A brightening of the Seyfert nucleus by a factor 10^ should 
therefore be explained by a sudden increase in the accretion 
rate. Tidal disruption s of objects by a mas s ive black ho l e (see 
Komossaetai]|2004t iBurrowset all 120111: ICenko et all 120121: 



Saxton et al.l2012h have been modelled and simula t ed by various 



aufliors (e.g.lReeJl988HE"vans & Kochaneklll989l:luimeJll999t 



Alexaride r & Kumai^ '2001'; [Alexander & Livid 120011: iLi et all 
2002t [Lpdato et al. 2009). The induced emission follows a 
power law decline with a characteristic slope of -5/3, corre- 
sponding to our observations. 

The peak of the observed 17.3-80 keV luminosity is 1.5 x 
10"*" erg s"'. Most of luminosity of a tidal disruption event is 
expected to be released in the soft X-rays (and absorbed by the 
torus in NGC 4845). We can evaluate the maximum soft X-ray 
emission allowed by the data by adding a soft thermal compo- 



Fig.8. The light curve of IGR J12580+0134 observed in 
the 17.3-80 keV energy band. Red squared points refer to 
INTEGRAL data, blue crosses show Swift and XMM-Newton ob- 
servations. The solid line show a model of the form (f - toY^^^, 
as expected for fallback of material after a tidal disruption event. 
The long-short dash and dotted lines indicate the predictions of 
simulations for the disruption of a sub-stellar object (y - 5/3; 
13 - 0.7) and, respectively of a star {y - 4/3; /3 = 0.65) 
(iGuillochon & Ramirez-Ruizll2013l) . 



nent to the model and determine for which flux such a com- 
ponent would have been detected. Figure [T] shows the increase 
in the x^ obtained depending on the soft component luminos- 
ity. For a temperature of kT - 0. 2 keV, which is expected in 
our case (see eq. 9 in IUlmeiJ[l999h . the soft component cannot 
be more than 10 times brighter than the hard X-ray component. 
The soft component could of course be brighter if it would peak 
at lower energies, which is however unexpected, especially for a 
low-mass black hole. 

Assuming that the thermal emission is ten times brighter than 
the hard X-ray emission, the tidal event luminosity reached the 
Eddington luminosity at maximum wit h Lflai-e /^Edd ~ 0.6, as ob- 
served in other tidal disruption events (iLi et al.ll2002l) . The total 
energy radiated is then of the order of 10^° ergs, which corre- 
sponds to the energy released by the accretion of 0.5 Jupiter mass 
(Mj). 

We fitted the decline of the IGR J12580+0134 X-ray Hght 
curve with the (f - t^y^^^ law where fp is the time of the initial 
tidal disruption (Fig. [8]l. The peak luminosity of the flare oc- 
curred on January 22, 201 1 according to MAXI and INTEGRAL. 
It turns out from our fit that the beginning of the tidal disruption 
occurred ^ 60-100 days before the peak of the X-ray flare, mark- 
ing the heating of the debris at t he vicinity of t he black hole. A 
similar delay was mentioned bv lLi et all ( |2002|) as the time Afi, 
since the most bounded material returns to the pericentre. 
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4.2. Comparison with tidal disruption simulations 

Detailed hyd rodynamic simulations of tidal disr uptions were 
performed by iGuillochon & Ramirez-Ruia (1201 3l) and their re- 
sults parametrised. The rate of mass falling on the black hole, 
M{t), the time of the peak accretion, fpeak, as well the decay 
power-law index, F, depend on a few parameters: the structure 
and mass of the disrupted object and the minimum distance to 
the black hole. 

We compared the simulations to our data assuming that 
the disrupted object is either a star or a sub-stellar ob- 
ject with polytropic index y of 4/3 or 5/3, respectively. 
We used the parametrisation incl uded into the Appendix of 
IGuillochon & Ramirez-Ruia (1201 3h for different impact param- 
eters jS = rx/rp, where rj, and rp are the tidal and the pericentric 
radii, respectively. We also assumed M*- /?* relations for the 
disrupted object valid for sub-stellar objects or stars according 
to lChabiier & Baraffe (2000) . 

The observational constraints, fpeak - 0.2 yr and the peak 
accretion rate Mpeak - 2.5 Mj/yr (under the assumption of a 
hard X-ray radiation efficiency of 10%), are sufficient to con- 
strain the mass of the disrupted object. For a black-hole mass of 
2.3 X 1O^M0, the mass of the disrupted object turns out to be 
either 14-16 Mj (fory = 5/3 and/3 = 0.6 - 1.9) or 10 - 15Mo 
(for y = 4/3 and /3 ^ 0.65). For a black-hole mass increasing 
to 1O*M0, the mass of the disrupted object would be 25 - 28Mj 
(up to 75Mj for lO^M©) for a sub-stellar object or 1 - 3M0 for 
a star (J3 ^ 0.7). The /3 parameter is not constrained in the case 
of a sub-stellar object, while a very naiTow range is required to 
disrupt a star. The fraction of the object reaching the black hole 
would be very small in the case of a stellar disruption. 

The hydrodynamic simulations also predict the detailed vari- 
ability of the accreted mass flow with time. In the case of a star 
(y - 4/3), less centrally condensed than a sub-stellar object, the 
accreted mass decreases much faster after the peak. This is il- 
lustrated in figure |8] where the observed flux behaves similarly 
to the predictions for the disruption of a sub-stellar object and 
contrasts with those obtained for a star. 

The hydrodynamic simulations indicate, therefore, that the 
tidally disrupted object was probably a 14 - 30Mj sub-stellar ob- 
ject and that about 10% of its mass has been accreted on a black- 
hole weighting no more than IO^Mq. We note that a slightly 
different equation of state could lead to a mass of the disrupted 
object as low as several My and that specific simulations of tidal 
disruption of sub-stellar objects have not yet been done. 

The decline in the hard X-ray flux observed in Fig. |8] after 
200-500 days can be understood as debris falling back on the re- 
maining object core an d decreasing the emission during th e late 
evolution (see fig. 8b in lGuillochon & Ramirez-Ruizl2013h . This 
would indicate that the disruption was not total. Another expla- 
nation could be that the corona disappears faster than the debris, 
possibly indicating a change in the geometry of the accretion 
flow with time. 
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Fig. 9. The 19.1 ks XMM-Newton power spectrum (blue his- 
togram) versus simulated 10 000 powers (red histogram with er- 
ror bars). The power spectrum of the simulated light curves are 
power law shaped with a slope of -0.9. Icr error bars are shown. 
An excess in the observed power near the log(frequency)= -3 
may suggest a presence of the QPO feature. 
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4.3. Power spectrum 

The power spectrum of the X-ray light curve obtained by XMM- 
Newton at the flare maximum (see Fig. |9]i shows an excess 
in the range 0.0008-0.004 Hz. This excess could be similar 
to a quasi-periodic oscillation (QPO) in the accretion flow. 
The frequency range of the excess corresponds to the inner- 
most stable circular orbit for a black hole of a few xIO^Mq 
dRemillard & McCUntoci2006l) . 



Fig. 10. Probability multiplied by number of trials (Mnais = 
10000) against x^ for simulated single and broken power-law 
models. The solid black curve points out expected probability of 
a perfect model. The best fit by a single power law is obtained 
with a slope = -1.3. Broken power laws with the low-frequency 
slope - -0.5, the high-frequency slope between -1.1 and -1,3, 
and with broken frequency = 0.001 Hz provide a slightly better 
description of the observed power spectrum. 
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To investigate the significance of this excess, we performed 
Monte Carlo simulations. We assumed various shapes of the 
power spectrum o f the source and generated lO'OOO light curves 
for each of them (iTimmer & Koeniglll995h . The spread of the 
power spectra of all simulated light curves were then compared 
to the shap e of th e observed power spectrum, as outlined in 
IUttlevetal.l(l2002h . 

We tried single power-law power spectra with slopes be- 
tween -0.7 and -1.5 and broken power-law models with break 
frequency frozen to 10"^ Hz, which differ by the low- and high- 
frequency slopes. Figure [TOlshows the;^'^ distributions obtained 
for the best models, together with the expected distribution. The 
deviation between the simulated and expected distribution indi- 
cates that the excess can be obtained by chance at a probability 
not greater than 5%. We conclude that a longer XMM observa- 
tion would have been very useful and note that an instrument like 
the Large Observatory For X-ray Timing (LOFT) would have 
detected several million X-ray photons during one ISCO period 
for this tidal event, which would have allowed a probe of the 
geometry of the debris falling towards the black-hole horizon. 

4.4. Tidal events frequency 

The density of stars in the centre of a well developed galaxy 
is around lO^Mg/pc^ in a region of 100 pc in radius and is 
comparable to what is found in globular cluster (iMarx & Pfaul 
Il992h . These dense regions allow stellar enco unters to create 
new g ravitationally bound systems or ejections (iDavies & Benzl 
Il995h . Stars in the vicinity of a supermassive black hole may 
escape t o infinity after encounter because their orbits are de - 
flected ( lAlexander & Kumail 120011: [Alexander & Livid l2001h . 
Other p rocesses may als o play important roles (see brief re- 
view b V I Alexandeij|20 1 2h . The tidal event rate in bulges are es- 
timated to be in the range 10"* - 10 "* yr"' galaxy ' growing 
to 10 ~ in case of chao tic stellar motions in a non- axisymmetric 
gravitational potential (iMerritt & Poonll2004h . A reasonable rate 
for the stellar tidal disruption events is ~ 10~^ yr~^ galaxy ~' 
dPonlev et al.ll2002Hmng & Merrittll2004t iGezari et al.ll2009h . 

The discovery of a population of un bound Jupiter-mas s ob- 
jects has recently be en mentioned by ISumi et aLl ( 1201 lb and 
iDelorme et al.l(l2012h . The number of free-floating massive plan- 
ets, which are kicked away from young sta r-forming regions 
(e.g. lOuanz et al.l [20101) . has been estimated dSumi et al.ll2011h 
as one to three times higher than the number of main-sequence 
stars. Because of the different equations of state, the tidal radii 
are similar for objects of IMq and IMj around a massive black 
hole. Tidal disruptions of planets could therefore be as frequent 
as star disruption. Planets could increase the rate of tidal events 
in galaxies a few times, but will generate lower luminosity events 
on average. 

We estimate that INTEGRAL could detect tidal disrup- 
tion flares with fluxes ten times lower than observed in 
IGRJ12580+0134 i.e. from galaxies up to a distance of 50 Mpc, 
for similarly faint events. The number of galaxies within D < 50 
Mpc is ~ 5000. The event rate detectable at hard X-rays (by 
INTEGRAL or Swift) could therefore be expected to be a few 
events every ten years. The INTEGRAL and Swift archives can 
be used to test these predictions. 
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